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ABSTRACT 

Within the framework of a two-fluid description possible pathways for the 
generation of fast flows (dynamical as well as steady) in the lower solar atmo- 
sphere is established. It is shown that a primary plasma flow (locally sub- 
Alfvenic) is accelerated when interacting with emerging/ambient arcade-like 
closed field structures. The acceleration implies a conversion of thermal and 
field energies to kinetic energy of the flow. The time-scale for creating reason- 
ably fast flows (> lOOkm/s) is dictated by the initial ion skin depth while the 
amplification of the flow depends on local j3. It is shown, for the first time, that 
distances over which the flows become "fast" are ~ 0.01 R s from the interaction 
surface; later the fast flow localizes (with dimensions < 0.05 Rs) in the upper 
central region of the original arcade. For fixed initial temperature the final speed 
(> 500 km/ s) of the accelerated flow, and the modification of the field struc- 
ture are independent of the time-duration (life-time) of the initial flow. In the 
presence of dissipation, these flows are likely to play a fundamental role in the 
heating of the finely structured Solar atmosphere. 
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1. Introduction 

In astrophysics (particularly in the physics of the solar atmosphere), plasma "flow" could 
be assigned at least two connotations: 1) The flow is a primary object whose dynamics bears 
critically on the phenomena under investigation. The problems of the formation and the orig- 
inal heating of the coronal structure, the creation of channels for particle escape, for instance, 
fall in this category, 2) The flow is a secondary feature of the system, possibly created as a by- 
product and/or used to drive or suppress an instability. Since the generation of flows which 
will eventually create the coronal structures (Mahajan et al. 1999; Mahajan et al. 2001) is 
the theme of this effort, the flows here are fundamental. 

By exploiting a simple two-fluid model in the solar context, several recent studies 
(Mahajan k Yoshida 1998; Yoshida k Ohsaki k Mahajan 2004) have revealed the breadth 
of phenomena made possible by the combined action of the flow-velocity and the magnetic 
fields. The flow-based approach will prove, perhaps, crucial in the study of solar corona, ob- 
servationally found to be a highly dynamic arena replete with multiple-scale spatiotemporal 
structures (Aschwanden et al. 2001a); the approach gains immense credibility with the dis- 
covery that strong flows are found everywhere - in the subcoronal (chromosphere) as well as 
in the coronal regions (see e.g. (Schrijver et al. 1999; Winebarger k DeLuca k Golub 2001; 
Aschwanden et al. 2001a; Aschwanden 2001b; Seaton et al. 2001; Winebarger et al. 2002; 
Feldman k Landi k Curdt 2003; Brynildsen et al. 2004) and references therein). Recent 
phenomenology strongly emphasizes that the solar atmosphere is an extremely inhomoge- 
neous (in all parameters) area in which small- and large-scale closed magnetic field struc- 
tures with different temperatures co-exist in nearby regions. For example, two-temperature 
coronal models constructed from SOHO/EIT observations indicate complicated magnetic 
topology and fine-scale structuring of corona (including Coronal Holes) (Zhang et al. 1999; 
Chertok et al. 2002). It is also clear that the mechanisms for energy transport and chan- 
neling of particles in Solar atmosphere are deeply connected with the challenging and ex- 
citing problems of the solar coronal heating and of the origin of the solar wind (SW) 
(Woo k Habbal k Feldman 2004). 

If flows are to play an important and essential role in determining the dynamics and 
structure of the solar corona, we must immediately face the problem of finding sources and 
mechanisms for the creation of these flows. Catastrophic models of flow production in which 
the magnetic energy is suddenly converted into bulk kinetic energy (and thermal energy) 
are rather well-known; various forms of magnetic reconnection (flares, micro and nano- 
flares) schemes permeate the literature (E.g. (Wilhelm 2001; Christopoilou, Georgakilas and 
Koutchmy 2001) for chromosphere up-flow generations). A few other mechanisms of this 
genre also exist: Uchida et al. (2001) proposed that the major part of the supply of energy 
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and mass to the active regions of the corona may come from a dynamical leakage of magnetic 
twists produced in the sub-photospheric convection layer; Ohsaki et al. (2001, 2002) have 
shown how a slowly evolving closed structure (modelled as a double-Beltrami two-fluid 
equilibrium) may experience, under appropriate conditions, a sudden loss of equilibrium 
with the initial magnetic energy appearing as the mass flow energy. Another mechanism, 
based on loop interactions and fragmentations and explaining the formation of loop threads, 
was given in Sakai and Furusawa (2002); the suggestion based on cascade of shock wave 
interactions was made in (Ryutova & Tarbell 2003). A more quasi-static mechanism for 
flow generation in sub-coronal regions taking into account the density in-homogeneity of 
the structures was given in (Mahajan et al. 2002). 

Before we embark on delineating the flow-generation mechanisms, we present additional 
evidence/speculation on their existence as well as their possible role in the processes taking 
place in the solar vicinity: 

1) Goodman (2001) has shown that the mechanism which transports mechanical energy 
from the convection zone to the chromosphere (to sustain its heating rate) could also supply 
the energy needed to heat the corona, and accelerate the SW. The coronal heating problem, 
hence, is shifted to the problem of the dynamic energization of the chromosphere. In the lat- 
ter process the role of flows is found to be critical as warranted by the following observations 
made in soft X-rays and extreme ultraviolet (EUV) wavelengths, and recent findings from 
the Transition Region and Coronal Explorer (TRACE): the over-density of coronal loops, 
the chromospheric up-flows of heated plasma, and the localization of the heating function 
in the lower corona (Schrijver, et al. 1999; Aschwanden et al. 2001a; Aschwanden 2001b). 

2) The connection/coupling of transient events like up-flows and different type jet- 
like structures to the photosphere dynamics was reported in numerous studies (see e.g. 
(Ryutova & Tarbell 2003) and references in). In (Liu et al. 2003) it was demonstrated that 
the eruptions of coronal mass ejection is triggered from the low solar atmosphere (photo- 
sphere/chromosphere) as seen in TRACE 1600A images and with SOHO Michelson Doppler 
Imager. The data of this latest research favor the idea that a catastrophic loss of MHD equi- 
librium can be the primary driving mechanism for the rapid ejection that has 3 important 
stages - a relatively stable equilibrium, a loss (fast, impulsive) stage, and the final rapid 
eruption (associated with substantive changes in the photospheric magnetic flux and white- 
light morphology). The results of (Socas-Navarro & Martinez Pillet & Lites 2004) suggest 
the coupling between magnetic fields and convective processes that pervades the solar pho- 
tosphere. The correlation between photospheric shear flows and flares is also reported in 
(Yang et al. 2004); several current models suggest that the former can be responsible for the 
energy build up of the flares. 
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3) In (Nitta & Oliver & Tylka 2003), the authors report on the low coronal signatures of 
major solar energetic particle (SEP) events focusing on flare-associated motions (observed in 
soft X-rays). It was underlined that these motions may provide an important link between 
small-scale energy release and large-scale explosive events; the existence of a continuum 
of acceleration timescales was pointed out. In (Magara & Longscope 2003) the detailed 
investigation of the dynamical behavior of emerging magnetic flux using three-dimensional 
MHD numerical simulation was carried out and it was shown that the emergence generates 
not only vertical but also horizontal flows in the photosphere, both of which contribute to 
the injection of the magnetic energy and helicity. The contributions of vertical flows are 
dominant at the early phase of flux emergence, while horizontal flows become a dominant 
contributor later. In (Falconer et al. 2003) it was shown that solar corona is mainly heated 
by the magnetic activity in the edges of the network flux clumps that are observed to be 
riddled with the fine-scale explosive events. They present that: (1) at the edges of the 
network flow clumps there are many transient sheared-core bipoles of the size and lifetime 
of granules and having transverse field strengths greater than ~ 100 G, (2) ~ 30 of these 
bipoles are present per supergranule, and (3) most spicules are produced by explosions of 
these bipoles. 

4) Recent observations also suggest that the energy for coronal heating is very likely 
a by-product of the outflow of heat from Sun's interior through the convection zone - the 
convection zone acts as a heat engine, converting some of the thermal energy into mechanical 
and magnetic energy, some of which enters the corona and dissipates into heat. There 
are only two obvious energy sources that could power significant flow generation in the 
chromosphere: the magnetic field (both large scale and short-scale including turbulence), 
and the thermal pressure of the plasma. We have already mentioned a few examples of the 
magnetically driven transient, but sudden flow-generation. A more quiescent pathway was 
studied in (Mahajan et al. 2002) showing the possibility of magneto-fluid rearrangement of 
a relatively constant kinetic energy (going from an initial high-density-low-velocity state to a 
low-density-high- velocity stage). The mechanisms based on the wave-energy transformation 
and instabilities can be operative at later stages of the flow evolution; these mechanisms could 
have additional importance for acceleration (Poedts & Rogava & Mahajan 1998). 

The main message then, is that to solve the coronal heating problem, the inclusion 
of processes taking place in the chromosphere and the transition region may be essential. 
In particular, one must take into account the different time-scale dynamical stages of the 
evolution of the primary flow as it passes through specific regions of solar atmosphere areas 
nested by varying scale ambient magnetic field structures. The dynamics of the flow must 
be understood. 
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In present paper we will show the possibility of flow acceleration/generation in the Solar 
atmosphere based on the dynamical two-fluid model suggested in (Mahajan et al. 2001). We 
will show that there exists an extremely fast stage (right at the lower chromosphere heights) 
giving rise to significant flow acceleration/generation; it is followed by a quasi-static stage in 
which the created fast flows are further accelerated via the magnetofluid coupling (depending 
on the region of the atmosphere the density could be constant or spatially varying). The 
detailed nature of accelerated flows will depend on the initial and boundary conditions. 



The physical model for flow generation/acceleration is a simplified two-fluid model. 
The plasma is quasi-neutral - - electron and proton number densities are nearly equal: 
n e ~ rii = n (V • j = 0), but the electron and the proton flow velocities are allowed 
to be different. Neglecting electron inertia, these are Vi and V e = (V — j/en), respec- 
tively. We assign equal temperatures to the electron and the protons so that the ki- 
netic pressure p is given by: p = pi + p e ~ 2nT, T = Tj ~ T e . The analysis can 
be readily extended later to the more realistic case of different temperatures for different 
species (McKenzie & Sukhorukova & Axword 1998). We understand that when solving the 
solar wind problem one should use the multi-fluid, multi-dimensional descriptions (see e.g. 
(Tu & Marsch 1997; Tu & Marsch 2001; Hollweg 1999) and references therein) but we how- 
ever, believe, that essential features of the primary flow-based physics of its acceleration can 
be captured with our basic model. Very near the photospheric surface, the influence of neu- 
trals and ionization (and processes of flux emergence etc.) would not permit just two-fluid 
approach. A little farther distance (Ar > 500 km) from the surface, however, we expect that 
there exist fully ionized and magnetized plasma structures such that the dynamical two-fluid 
model will capture the essential physics of flow generation. 

The dimensionless two-fluid equations describing the flow-field interaction processes 
can be read from (Mahajan et al. 1999, 2001): 



2. 



Model 



d_ 

dt 






V-b = 0, 

d 

— n + V-nV = 0, 



(3) 



(4) 
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^(2T)+V(q, + q e ) 



+ ^« (V x b) • VT - — (V x b)V(nT) + £^ - 



(5) 



where the notation is standard with the following normalizations: the density n to n at 
some appropriate distance from the solar surface, the magnetic field to the ambient field 
strength at the same distance, and velocities to the Alfven velocity Vao- The parameters 
r A0 = GMq/V% R q = 2/3o/r c0 , «o = Ko/Rq, Po = c 2 s0 /Vjy are defined with n , T , B . 
Here CgQ IS db sound speed, R Q is the solar radius, r c0 = GM Q /2c 2 R Q , A^o = c/ooiq is the 
collisionless skin depth, fi(n, T) is ion kinematic viscosity and q e and g« are electron and ion 
dimensionless heat flux densities, Eh is the local mechanical heating function and Er is the 
total radiative loss. We note that the full viscosity tensor relevant to a magnetized plasma 
is rather cumbersome, and we do not display it here. Just to have a feel for the importance 
of spatial variation in viscous dissipation, we display its relatively simple symmetric form. 
It is to be clearly understood that this version is meant only for theoretical elucidation 
and not for detailed simulation. We also note that in general, Hall current contributions 
are expected to become significant when the dimensionless Hall coefficient satisfies the 
condition: ao > rj, where r) is the inverse Lundquist number for the plasma. For a typical 
coronal plasma as well as for low chromosphere and transition region (TR) this condition is 
easily satisfied (a is in the range 1CT 10 — 10~ 7 for densities within (10 14 — 10 8 ) cm~ 3 and 
rj = c 2 / (AkVaoRqv) ~ 1CT 14 , where a is the plasma conductivity). 

To establish the relevant parameter regime, we resort to recent observational data 
(e.g. (Goodman 2001; Aschwanden et al. 2001a; Socas-Navarro & Sanchez Almeida 2002) 
and references therein). At ~ (500 — 5000) km, the observations yield the average plasma 
density and temperature to be respectively n ~ (10 14 — 10 11 )cm~ 3 , and T ~ (1 — 6) eV. 
For simplicity, we have assumed T e = Tj = T. The information about the magnetic field is 
hard to extract due to the low sensitivity and lack of high spatial resolution of the measure- 
ments coupled with the inhomogeneity and co-existence of small- and large-scale structures 
with different temperatures, (observational evidence of small scale mixtures of weak and 
strong fields (Socas-Navarro & Lites 2004)) in nearby regions. The observation of pixel-to- 
pixel variations in the magnetic field indicates that small-scale (sub-pixel) distribution of 
fields changes considerably at larger spatial scales (Socas-Navarro & Sanchez Almeida 2003; 
Socas-Navarro 2004). At these distances we have different values for the network and for 
the internetwork fields: (i) The network plasmas have typically short-scale fields in the range 
Bq ~ (700 — 1500) G, have more or less n ~ const . (ii) The internetwork fields are generally 
weaker (with some exceptions (Socas-Navarro & Sanchez Almeida 2002)) — B Q < 500 G, 
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and are embedded in larger-scale plasma structures with n ^ const . For different classes of 
magnetic field structures different scenarios may be operative. 

3. Acceleration of particle flows analysis of the stages 

In our investigation we shall assume that the processes that generate the primary flows 
and the primary solar magnetic fields are independent. The plasma flows begin to interact 
with the ambient field at time t=0. The choice of initial conditions for our numerical work 
is guided by the observational evidence presented in the introduction. Our approach is 
consistent with that of Woo, Habbal and Feldman (2004) who have argued that the flow of 
the solar wind is influenced by the closed field structures stressing the self-consistent process 
of acceleration and trapping/heating of plasma particles in the finely structured atmosphere. 
We will dwell, in this paper, on the representative problem of the trapping and acceleration of 
the primary flow impinging on a single closed-line structure. The simulation was performed 
for a variety of initial and boundary conditions and essential aspects of the typical results 
will be presented below. 

3.1. Dynamical generation of fast flows 

The general set of Eqs. (1-5) was solved numerically in Cartesian Geometry for 2.5 
Dimensions (d y = 0). Note that the 2.5D Cartesian nature of our code does not allow us to 
explore large distances from the surface due to interference with the boundaries. Fortunately 
that does not translate into a serious shortcoming because much of the action is found to 
be limited to regions very close to the surface; the simulation results, therefore, are quite 
trustworthy in the revelation of the basic processes of interest. In carrying out the simulations 
an important assumption was made: the diffusion time of magnetic field is longer than the 
duration of the interaction process. 

A short summary of our numerical methods is in order. We use the 2.5D version of 
Lax-Wendroff finite difference numerical scheme along with applying the Flux-Corrected- 
Transport procedure (Richtmyer & Morton 1967; Zalesak 1979). The predictor-corrector 
type of approximation was used. Equation (3) was replaced by its equivalent for the y- 
component of the vector potential to ensure the divergence-free property of the magnetic 
field. The equation of heat conduction was treated separately by the alternate direction 
implicit method with iterations. Transport coefficients for heat conduction and viscosity 
are taken from Braginski, 1965. In the code, the Bremsstrahlung radiation accounts for E R 
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(Mahajan et al. 2001), using a somewhat modified formula assuming it to be 2 times greater, 
E R = 2 • E Br = 2 • 1.69 ■ 10~ 25 • n 2 ■ T 1/2 ■ Z 3 erg/cm 3 s, (Z — 1). Since we were exploring a 
particular heating mechanism suggested there, no external heating source Eh was needed. 
A numerical mesh of 280 x 220 points was used for computation. The corresponding scheme 
is characterized by second order accuracy with respect to the chosen grid. 

The latest observations support the idea that the coronal material is injected discontin- 
uously (in pulses or bunches, for example) from lower altitudes into the regions of interest 
(e.g., spicules, jet-like structures). A realistic simulation, then, requires a study of the inter- 
action of both temporally and spatially nonuniform initial flows with arcade-like magnetic 
field structures. These "close to the actual" cases represent more vividly the dynamics of 
the flow acceleration process. Below we study the dynamics of spatially-temporally localized 
flow (initially a Gaussian, Fig.-s 1,3) entering the region nested with varying scale arcade- 
like closed field line structures. For better visualization of the results we take the symmetric 
case. The flow is assumed to be initially weak (|V| 0maa , <C C s q). The initial ambient 
magnetic field was modelled as a single 2D arcade with circular field lines in the x-z plane 
(Fig. 2 for the vector potential/flux function). The arcade field attains its maximum value 
B max (x ,z — 0) = B 0z at x at its center, and is a decreasing function of the height z 
(radial direction). This field was assumed to be initially uniform in time. When doing so, we 
choose the parameters to satisfy the observational constraint that, over a period of some tens 
of minutes, the location of the trapping/acceleration must have a relatively smooth evolu- 
tion. The final shape and location of the structure of the associated B(r, t), for example will 
be naturally defined by its material source, by the process dynamics, and by the initial field 
B (r, t). We use the following representation for the magnetic field: B = VA y + B z z and 
for given geometry A(0; A y ; 0); b = T5/Bq z ; b x (t, x, z ^ 0) ^ 0. From numerous runs on the 
flow-field evolution, we have chosen to display pictorially the results corresponding to the 
following initial and boundary conditions: B 0z = 100 G and flow parameters: V max (x , z = 
0) = V 0z = 2.18 • 10 5 cm/s; n 0max = 10 12 cm" 3 ; T(x,z = 0) = const = T = 10 eV. The 
background plasma density n\> g = 0.2n 0max . In simulations n(x,z,t = 0) = n/no max is an 
exponentially decreasing function of z. Experienced gained after numerous runs, revealing 
that the processes under study are localized within a small area of interaction, we settled on 
the following boundary condition, d x JC(x = ±oo, z,t) = which was used with sufficiently 
high accuracy for all parameters /C(A, T, V, B, n) . Guided by observations we assume that 
the initial velocity field has a pulse-like distribution (Fig. 3) with a time duration (life-time) 
to > 50 s. 

To illustrate the acceleration of initial flows (extremely weak), we have modelled several 
cases with different initial and boundary conditions. The dynamical picture is strongly de- 
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pendent on the relative strengths of the initial flow pressure and the magnetic field strength. 

Our typical representative example is the evolution of a symmetric weak up-flow with 
its peak located in the central region of a single closed magnetic field structure (location of 
field maximum B 0z = 100 G) (Figs. 1-3). Figs. (4-8), in which we give the x — z contour 
plots of all the relevant fields (A y ; |b|; n; |V|; T), contain the essence of the simulation. 
We find that the acceleration is significant in the vicinity of the field-maximum with strong 
deformation of field lines and energy re-distribution. In this very region, the simulations 
show cooling of the flow with serious density redistribution: part of the flow is trapped in 
the maximum field localization area, accumulated, cooled and accelerated. The accelerated 
flow reaches > lOOkm/s value in less than 100 s (in agreement with recent observations 
(Ryutova & Tarbell 2003; Seaton et al. 2001; Schrijver et al. 1999) and references therein). 
The accelerated flow is decoupled from the mother flow, and is localized in a distinguishable 
region with dimensions < 0.05 R s starting at a distance ~ 0.01 R s from the interaction 
surface. The time for reaching the quasi-equilibrium parameters is determined by the initial 
and boundary conditions (this conclusion is general for all cases). 

In a stationary analysis to be presented in the next sub-section, we will attempt to 
derive the characteristic steady state parameters (like the distance from the surface) of the 
simulated system. 

Extensive simulation experiments show that, when viscosity and heat flux effects are 
included, the flow acceleration evolution parameters depend strongly on a , the parameter 
measuring the strength of the Hall term in the two-fluid equations. A very interesting and far- 
reaching result is that the final parameters of the accelerated flow are practically independent 
of the initial flow-characteristics (Fig.-s 4-8); only the initial fast stage of acceleration up to 
~ 200 km/ s is slightly different for different primary flows. Simulation results for 2 different 
initial life-times of the flow (£ = 1000 s - left panel and to — 100 s - right panel in Fig.-s 
4-8) illustrate this feature. 

We also found that at some critical time, the solutions split into two parts; all fields, 
the magnetic (Fig.-s 4,5), the density (Fig. 6), the velocity (Fig. 7) and temperature (Fig.8) 
exhibit bifurcation. This process persists for different initial conditions. In Fig. 9 we give 
time evolution plots of the maximum values of all fields (A y , |b|, b p , b z , n, |V|, V p , V z , T) 
for a pulse-like flow interacting with a single arcade-structure for different initial life-times 
(t ) of the flow (t = 100 s{black); t = 1000 s(red); t = 2000 s (blue); t -> oo (green)). In 
Fig. 10 the same maximum values of all fields are plotted versus the initial life-time (t ) of 
the flow for different time-frames (t = 200 s(black); 500 s(red); 1000 s{blue); 1500 s(green); 
2000 s(lightgreen); 3000 s (rose)). 
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These pictures clearly demonstrate that the accelerated mother flow bifurcates into 2 
separate fast daughter-flows (after an initial acceleration stage) modifying significantly the 
original arcade structure. The characteristic fields undergo rather similar dynamics for flow 
pulses with different initial life-times. It should be emphasized that now the maxima of these 
parameters are localized not along the initial B-maximum but on both sides of it and shifted 
along height (in the localization areas of each accelerated daughter flow with newly created 
B-maxima and currents). After the initial acceleration stage, the magnetic energy maxima 
remain practically unchanged up to some "blow-up" time (> 2000 s) at which the gradients 
become too steep and the simulation results cease to be meaningful. The same result holds 
for the maxima of the transverse and parallel magnetic field energies (with < 10% accuracy). 
For fixed To(/3q) the maximum values of each parameter (local in space) exhibit practically 
similar dynamics (independent of the initial flow life-time) reaching similar numbers at near- 
critical time. This picture persists for different initial T (/3o). Testing the conservation of the 
total energy of the system as it evolves in time also shows that the simulation results can be 
trusted only up to the blow-up time; as one approaches this time the energy conservation 
no longer holds. To study longer time dynamics, the code will need improvement. 

We will soon offer a possible explanation of these results through a simple equilibrium 
analysis. 

We are now in a position to list the most interesting and distinguishable new results 
found in a 2.5D simulation of the two-fluid equations (containing various dissipative and 
short-scale effects) solved for different initial and boundary conditions: 

1. A primary flow, even with a very slow initial speed (V 0z ~ 1 km/s locally sub-Alfvenic) 
is accelerated when it interacts with an arcade-like closed magnetic field structure. The 
effect is strong in the strong field region (initially the arcade center) . This is a common 
feature independent of the arcade-characteristics, and the shape of the initial flow. 

2. When viscosity and heat-flux are ignored, the time needed for the flow to acquire 
reasonable amount of energy is practically infinite. This is probably due to the fact 
that without dissipation, the energy transfer through the short-scales introduced by 
the two-fluid effects is not operative unless special conditions for catastrophic processes 
pertain. 

3. For realistic a (measuring the strength of the Hall term) the heat flux and viscosity 
effects cause a re-distribution of magnetic, flow kinetic and thermal energies in the 
arcade region in reasonable times ~ 100 s . 
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4. During the redistribution, the arcade field is modified; the thermal- and magnetic 
field- energies are converted to flow energy locally. The time-scale for generating a 
reasonably fast flow (V 0z > lOOkm/s) is dictated by a . For a given initial T (/3 ), 
the larger the , the faster the flow generation. The density is non-uniformly 
redistributed within the arcade span. 

5. At some specific "critical time" < 1000s the accelerated flow bifurcates into two 
separate fast flows. At this moment the arcade is also split in two, each with its share 
of the accumulated particles. Two fast spicule-like structures, carrying vorticity and 
current, are decoupled from the mother flow. Their densities are similar to the initial 
density of the mother-flow. 

6. Initially the amplification of the flow depends on the flow fio , the ratio of the thermal 
and the magnetic field energy. 

7. The distance from the interaction surface where the bifurcation occurs is ~ 0.01 R s 
. It is interesting to mention that this height is lower than the heights of the base 
of a typical hot coronal structure ((Schrijver et al. 1999; Mahajan et al. 2001) and 
references therein) and it seems to be comparable to the latest observational find- 
ings (Brynildsen et al. 2004; Woo & Habbal & Feldman 2004). Initially the fast flow 
localizes in the center of the original arcade, starting from this distance. After the bi- 
furcation several flows appear with distinguishable dimensions (< 0.05 R s ) practically 
on similar heights. 

8. For fixed initial To, the final speed of accelerated flow and the picture of the modified 
field structure are independent of the initial flow life-time. This result seems extremely 
important in connection with the observed flows in the lower atmosphere. At t > 1000 s 
velocities reach ~ 500km/s or even greater (< 800km/s) values. Such result persists 
for different T (/9o) • 

We note here that at any quasi-equilibrium stage of the acceleration process, the nascent 
intermittent flows will blend and interact with pre-existing varying scale closed field struc- 
tures (recall the fine structure of the solar atmosphere); the "new" flows could be trapped by 
other structures with strong/weak magnetic fields participating in creating different dynam- 
ical scenarios: heating of the new structure (Mahajan et al. 2001) could result, or an escape 
channel could be created (Mahajan et al. 2003; Woo & Habbal & Feldman 2004). Instabil- 
ities, generation of waves could also be triggered. 
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3.2. Quasi-equilibrium pathway of flow acceleration due to magneto-fluid 

coupling — restrictions and analysis 



Both the observational evidence and the results of dynamical simulation point out that 
a typical solar structure passes through a quasi-equilibrium stage (possibly even a series of 
quasi-equilibria) before it reaches the final explosive or distortion/deformation stage leading 
to particle escape. Let us try to understand the physics of these quasi-equilibrium structures 
in terms of equilibrium two-fluid equations. We analyze the simplest two-fluid equilibria 
with T = const — > n _1 Vp — > TV In n (generalization to a homentropic fluid: p = 



const • n 7 is straightforward and was performed in numerical work (Mahajan et al. 2002)). 
The dimensionless equations describing the model equilibrium can be written as: 



1, 



-V x b x b + V 



n 



— - A) In n - — ) V x (v x V) = 0. 



V x 



(v-%xb)xb]= 0, 

V • (nV) = 0, 
V-b = 0, 



(6) 

(7) 
(8) 
(9) 



where b = B/i?o and the following normalizations were used: n — > no - the density at 
some appropriate distance from the solar surface (> 2000 km), B — > B - the ambient field 
strength at the same distance, |V| — > Vao and the dimensionless parameters are defined with 
n , T , B taken at the same distance. In the non-dissipative limit, the system allows the 
well-known double Beltrami solutions : 



b + a V x V = d n V, 



a n 



V 



— V x b 

n 



(10) 



where a and d are the dimensionless constants related to ideal invariants: the magnetic 
hi — /(A • b) d 3 x and the generalized h 2 — /(A + V) • V x (A + V)d 3 x helicities 
(Mahajan & Yoshida 1998; Mahajan et al. 2001). Substituting (10) into (6)-(8) one obtains 
the Bernoulli Condition 

'2/3 r cn „ , V 2 ' 



- Pq In n - — 



0, 



(11) 



relating the density with the flow kinetic energy, and solar gravity. 



Equations (6), (10), (11) represent a close system. They may be easily manipulated to 
yield an alternative form (g(r) = r c0 /r) 



-^V xVxV + a Vx 
n 



d I n V 

an 



d -\ V 

a 



0, 



(12) 
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«nV X 




< bj + a V x ^— - dj b + ^1 - b = 
= exp(-[ 2flb -^--2^ + — j). 



(13) 



(14) 



Equations (12), and (14) can be solved for the density and the velocity field V and 
then b could be determined from (10). 

In the Solar atmosphere one observes quasi-equilibrium magnetic structures with both 
homogeneous (practically anywhere) and inhomogeneous (especially in the Chromosphere 
and TR) densities. By invoking appropriate variational principles, one can show that the 
generic double Beltrami class of equilibria are accessible in all cases of interest: constant 
density, constant temperature, or when the plasma obeys an equation of state. Maximum 
analytical headway, however, is possible for constant density. In that case the Bertrami- 
Bernoulli system consists of a set of linear equations and has two well-defined scales of 
variation. Non-constant density does not lead to a linear chain (see (12), and (14)), but 
allows phenomena peculiar to nonlinear systems. It is the latter class of systems that we 
will deal with now. 

We will now calculate the amplification conditions for inhomogeneous density flows 
in the chromosphere. We restrict to a one-D variation (along the height Z) and choose 
the constants a ~ d = 100 so that (a — d)/ad ~ 10~ 6 . This choice insures that two 
homogeneous Beltrami scales will be vastly different. Detailed algebraic derivation of the 
approximate formulas used below can be found in Appendix 1. 

The principal results of Appendix 1 are that if n^> (ad)~ l (density fall in the region 
of interest is not more than 3 orders of magnitude), then 

1) the transverse components of magnetic field vary keeping b 2 + b 2 = 6qj_ = const. 

2) The density and the velocity fields are related approximately by | V" | 2 = l/d 2 n 2 so 
that the magnetic energy does not change much, |b| 2 = const to leading order. 

3) The Bernoulli condition transforms to the defining equation for density: 



We notice that for the density to drop with height, it has to be larger than n min = (2(3 ) 1 ^ 2 d. 
The existence of n min also implies via V 2 = l/d 2 n 2 that the maximum allowed velocity is 




(15) 



1 



(2A,) 1/2 . 



(16) 



max 



dn 



■mm 
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As one approaches the singularity at n — n min , the spatial variation of density (and in 
particular of the velocity) becomes very large. In such a region of the steep fall in density 
and rise in velocity, the time-independent dissipationless approach will not be valid. The 
Bernoulli equation (15), however, clearly reveals the origin of the very fast first stage of 
dynamical acceleration found in the simulations. From Eq.(15) we also see that the distance 
over which the catastrophe appears is determined by the strength of gravity, g(z) . Eventual 
amplification of the flow is determined by the local value of (3 . These simple consequences 
of the Bernoulli equation explain one of the most important findings of the simulation: for a 
fixed initial temperature, the final characteristic parameters of the accelerated flow (quasi- 
equilibrium after the fast stage of acceleration) do not depend on its initial state. For these 
gross features of the system, the value of as long as it is finite, is also quite irrelevant, it 
just determines how fast the transverse components of magnetic and velocity fields oscillate. 
However when dissipation is present, through the mediation of short-scale physics, plays 
a crucial role in the acceleration/heating processes. 

In connection with this result it is interesting to mention that according to latest obser- 
vations on the quasi-equilibrium coronal loops, the so called quasi-periodic intensity oscilla- 
tions are found to propagate upwards with speeds of the order of the (adiabatic/isothermal) 
coronal sound speed ((De Moortel & Parnell $ Hood 2003) and the references therein). 

For structures with (n = const), there are two distinct scenarios for eruptive events in 
the current framework : (1) when a "slowly" evolving structure finds itself in a state of no 
equilibrium and (2) when the process of creating a long-lived hot structure is prematurely 
aborted; the flow shrinks/distorts the structure that suddenly shines and/or releases energy 
or ejects particles. The latter mechanism requires a detailed time-dependent treatment. 
The semi-equilibrium, collisionless magnetofluid treatment pertains only to the former case 
(Ohsaki et al. 2001; Ohsaki et al. 2002). In the references cited, the conditions for catas- 
trophic transformations of an original DB (double Beltrami state) were investigated. It was 
shown that when the total energy of the original state exceeds a critical value, the DB equi- 
librium suddenly relaxes to a single Beltrami state corresponding to the larger macroscopic 
scale; at the transition, much of the magnetic energy |b| 2 of the original state is converted 
to heat /flow kinetic energy. 

4. Summary of the results and Conclusions 

We have developed a 2.5 Dimensional dynamical code for two-fluid equations. The two 
fluid equations contain the Hall term (a ^ 0), the ion vorticity, heatflux and viscosity 
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effects. We have used the code for a systematic study of particle acceleration and energy re- 
distribution phenomena associated with the interaction of a primary plasma flow with closed 
field-line magnetic structures. We also developed simple analytical arguments to explain and 
understand essential features of the simulation results. The simulation and analytical effort 
have led us to several far-reaching results for the understanding of the solar atmosphere. 
Even at the cost of some repetition, we list the most important ones : 

(1) A primary plasma flow (locally sub-Alfvenic) is accelerated when it impinges on an 
emerging/ambient arcade-like closed magnetic field structure. The effect is strong in the 
strong field region. It is found that the final state of the flow is quite insensitive to the 
details of initial and boundary conditions; the latter simply dictate the time-scale at which 
significant flow-energy is generated. 

(2) It is shown that there is a redistribution of magnetic, flow- kinetic and thermal energies 
in the arcade region so that the original arcade field is modified, and thermal and field en- 
ergies are converted to flow energy. The time-scale of the fast (> lOOkm/s) flow generation 
is dictated by ao, the measure of the Hall term . 

(3) It is found that at some specific time the accelerated flow bifurcates into 2 separate fast 
flows with an accompanying split of the arcade each containing its share of the accumulated 
particles. 

(4) Initially the amplification of the flow depends on /3q as proven by the ID analysis; it is 
shown for the first time that the distance on which it happens is ~ 0.01 R s (independent of 
ao) from the interaction surface. Later the fast flow localizes (with dimensions < 0.05 Rs) 
in the upper center of the original arcade. 

(5) It is shown that for fixed initial T the final speed (> 500 km/ s) of the accelerated flow, 
and the shape of the modified field structure are independent of the initial flow life-time. 
Many of these parameters can be approximately calculated by analysis. 

We have shown possible pathways for both dynamical and steady generation of fast 
flows. The cold flows originating, for example, in the lower chromosphere acquire energy as 
they meet and interact with emerging/ ambient magnetic fluxes; the trapping of an ionized 
> 3 eV plasma by network/ inter-network structures takes place at the same time. In the 
presence of dissipation, these flows are likely to play a fundamental role in the heating of 
the finely structured solar atmosphere. The explicit purpose of this paper, however, was to 
demonstrate the generation of flows in the lower atmosphere feeding on the ambient magnetic 
energy. The flows, in turn, provide a steady and assured source of matter and energy for the 
formation and primary heating of the corona as well as for the creation of the solar wind. The 
agreement of our preliminary results with the observation data lends credence and promise 
to attempts, based on the exploitation of sub-coronal flows, to tackle unresolved problems 
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like the coronal heating and origin of the solar wind. We believe that the chromospheric mass 
outflows, spicules, explosive events in chromosphere, micro- and nano-flares, large coronal 
flares, erupting prominences and CMEs may happen separately but can also be parts of a 
more global dynamic process of coronal specific regional formation. 
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Appendix 1 — Equilibrium analysis of particle acceleration for non— uniform 
density case due to Magneto-fluid coupling 

Let's rewrite DB equations (10) in following way: 

a V x b = --b + n V, a V x V = -b + dnV , (1) 

a 



Let's define a vector: 

then (1) may be written as: 
where 



Q=[ y )' (2) 



a V x Q = M Q, (3) 



(4) 



M can be diagonalized by a similarity transformation: 

where [A 2 — (dn — a -1 ) A + n (1 - d cr 1 ) = 0] A± = ±[(dra - a" 1 ) ± + a" 1 ) 2 - 4n ] 

are standard roots. S is found to be (n is a slowly varying parameter, see the Bernoulli 
condition - V 2 and g are slowly varying): 



1, -(A+H-a- 1 ) 

1, -(X-+CL- 1 ) 



(6) 



Then, if density fall is at a much slower rate than the slow scale of the Beltrami system 
\_/a , ), rewriting (3) as: 



a V xSQ= (SMS-^SQ = ^ ^ + ' A ° ^ 5Q , 



(7) 



one finds: 



each obeying its own independent (fully de-coupled) equation: 



VxQ± = ^g±. (9) 

a 
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Let's find the amplification conditions for flows (say in the chromosphere, where a ~ d = 100 
so that (a — d)/ad ~ 10~ 6 ). Assuming (this is found to be a restriction) n ^> (ad)^ 1 
density fall is not more than 3 orders of magnitude, then 

A + ~ dn , A_ ~ — . (10) 

a d 

Notice, that for realistic solar atmosphere parameters (chromosphere, TR, corona) a ~ 
10~ 9 — 10~ n and the fast Beltrami scale \+/a ~ 10 11 — 10 13 is very oscillatory and its 
amplitude must go to zero. This gives a relation between the velocity and the magnetic field; 



Q + = b-(dn-a~ 1 )V ~b-dnV = 0, (11) 
and the approximate equation for the pertinent solution takes the form 

V x Q_ = Q_ with Q_ = b - ^ ~ b. (12) 

ad at) d 

Let's consider a ID problem (Z along height, b — 1 when normalized). Eq.(12) leads to: 

d 



dz 



(bl + b 2 y )=0 =► bl + bl = bl ± . (13) 



Then, using eq.(ll), one has: V% + Vy = b^ ± /d 2 n 2 . From Continuity Equation and DB 
condition: V z = Vo z /n ~ bo z /dn . Thus, 

V " = 1^-2- ( 14 ) 

Eq.(14) converts the Bernoulli condition (T = const) to: 

1 \ dn 



- 2 *" 2 + ^JS = " 39 - < 15 » 

Notice, that maximum allowed velocity for this mechanism is (compare with the condition 
(10) of (Mahajan et al. 2002)): 

\V max \ = j^— = (2(3 f 2 . (16) 
Analysis gives similar results for varying temperature (T = < /i < 1). 
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Fig. 1. — Initial symmetric profiles of the radial velocity V z , and density n. The respective 
maxima (at x—0) are ~ 2 km/ s and I0 12 cm -3 . 

Fig. 2. — Contour plots for the y- component of vector potential A (flux function) in the 
x — z plane for a typical ambient arcade-like solar magnetic field (initial distribution). The 
field has a maximum B max (x = 0, z = 0) = 100 G . 

Fig. 3. — The original pulse is limited in time. A time plot of V Zymax (t, z — 0) corresponding 
to the shape V z (t, z = 0) = V 0z sin(nt/to); V z (t > t ) = 0. The parameter t (1000 s for this 
pulse) can be interpreted as the "life-time" of the pulse. 

Fig. 4. — x — z contour plots at various time-frames: t = 200 s; 500 s; 1000 s; 1500 s; 
2000 s; 3000 s for the dynamical evolution of A y for flows with two different initial life- 
times. The spatially and temporally inhomogeneous (type displayed in Fig.l, Fig. 3) primary 
flows are accelerated as they make their way through the magnetic field with an arcade- 
like structure (Fig. 2). The realistic viscosity and heat-flux effects as well as the Hall term 
(a = 3.3 ■ 10~ 10 realistic) are included in the simulation. Left panel corresponds to the case 
of initial flow life-time: to = 1000 s, right panel — for t = 100 s. There is a critical time 
(< 1000 s) when the accelerated flow bifurcates in 2; the original arcade field is deformed 
correspondingly. 

Fig. 5. — x — z contour plots for the dynamical evolution of |b| exactly following the pattern 
of Fig.4. After the bifurcation (read caption of Fig.4), strong magnetic field localization 
areas, carrying currents, are created symmetrically about x — 0. 

Fig. 6. — x — z contour plots for dynamical evolution of density n exactly following the 
pattern of Fig.4. Post-bifurcation daughter flows are localized in the newly created magnetic 
field localization areas. The maximum density of each daughter flow is of the order of the 
density of the mother-flow. Daughter-flows have distinguishable dimensions ~ 0.05 R s 

Fig. 7. — x — z contour plots for the dynamical evolution of |V| exactly following the 
pattern of Fig.4. The initial flow, locally sub-Alfvenic, is accelerated reaching significant 
speeds (> lOOkm/s) in a very short time (> 100s). The effect is strong in the strong field 
region (center of the arcade). At t > 1000 s, the velocities reach ~ 500km/s or even greater 
(< 800km/s) values. The distance from surface where it happens is > 0.01 R s . 

Fig. 8. — x — z contour plots for dynamical evolution of temperature T exactly following 
the pattern of Fig.4. In the regions of localization of the daughter flow there is a significant 
cooling while the nearby regions are heated. 

Fig. 9. — Dynamical evolution of the characteristic fields (their maximum values) , |b|; b x ; b y ; 
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n; |V|; V x ; V y ; T, defining the interacting flow-magnetic field system (their x — 
z contour plots are shown in Fig.-s 4-8) for different initial flow life-times (to = 
100 s (black); 1000 s (red); 2000 s (blue);oo (green)). The code ceases to be dependable 
for times at which very steep gradients emerge; the blow-up time for this simulation is 
(< 3000 s) 

Fig. 10. — Maximum values of |b|; b x ; b y ; n; |V|; V x ; V y ; T (their x — z contour plots are 
shown in Figs. 4-8) versus the initial life-time (to) of the primary outflow for different 
time-frames (t = 200 s (black); 500 s (red); 1000 s (blue); 1500 s (green); 2000 s (lightgreen); 
3000 s (rose)). The code ceases to be dependable for times at which very steep gradients 
emerge; the blow-up time for this simulation is (< 3000 s) 
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